Fermi National Accelerator Laboratory FERMlLAB-Pub-02/028-T 

hep-ph/0202174 
February 1, 2008 

Increased Yield of tibb at Hadron Colliders in Low-Energy 

Supersymmetry 



(N 
O 

o 

(N 



o 

Oh 



Adam K. Leibovich and David Rainwater 
Theory Dept., Fermi National Accelerator Laboratory, Batavia, IL, USA 



Abstract 



Light bottom squarks and gluinos have been invoked to explain the b quark pair 
£SJ ■ production excess at the Tevatron. We investigate the associated production of tibb at 

hadron colliders in this scenario, and find that the rates for this process are enhanced 
. over the Standard Model prediction. If light gluinos exist, it may be possible to detect 

them at the Tevatron, and they could easily be observed at the Large Hadron Collider. 

o 

The bottom quark pair production cross section measured at the Fermilab Tevatron exceeds the 
theoretical prediction by about a factor of two The state of the art theoretical prediction is 
currently next-to-leading (NLO) order in QCD. While the NLO corrections are large, it is possible 



i \ that the measured excess is due to new physics 0. 

Berger et al. |§ propose a solution to this puzzle based on low-energy supersymmetry (LESS) [§J. 
i-Ch ■ In particular, they proposed the existence of a light gluino, g, with mass rrig ^12—16 GeV, which 
>• ■ decays to a bottom quark and light bottom squark, b\, of mass « 2 — 5.5 GeV. The bottom squark 
^ ! is either long lived or decays hadronically. It is argued that this scenario is not yet ruled out by other 
^ \ experiments |5|||. 

One way to test the above scenario is via observation of the final state tibb. With LESS, there is 
the new production channel tigg, with the gluinos decaying immediately to bb\. If the bottom squark 
decays hadronically, the decay products are typically merged with those of the associated b quark 
jet. If instead the bottom squarks are long-lived, the event signature is still two additional bottom 
quarks in top quark pair production. One would expect the rate for this new channel to be larger than 
the Standard Model (SM) tibb rate due to the large Casimir. Looking for LESS in this manner has 
many nice properties. Unlike for T |] or B § meson decays, we do not need to worry about non- 
perturbative physics. Furthermore, the scale dependence for top quark pair associated production is 
much smaller than for bottom quark pair production. Finally, unlike many supersymmetric searches, 
with this production channel there is almost no model dependence. As the production coupling 
involved is that of QCD, a s , to a first approximation the only LESS model parameter which enters 
at leading order is the gluino mass, m- g (inclusion of top squarks leads to additional diagrams, but 
these are suppressed primarily due to the additional heavy propagators). One could also choose to 
examine other associated production processes, such as Zgg. The larger scale uncertainty could be 
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compensated by the comparatively larger cross section, but the processes are not purely QCD, so we 
do not consider them here. 

We perform leading order, parton level Monte Carlo calculations of the SM ttbb and LESS tigg 
production cross sections, for gluino masses in the range m g = 12 — 16 GeV. Cross sections are 
calculated for both pp collisions relevant to Run II of the Fermilab Tevatron, y/s = 2.0 TeV, and for 
pp collisions at the CERN Large Hadron Collider (LHC), y/s = 14.0 TeV. Decay of the top quarks 
is included at the matrix element level, to determine the efficiency of realistic kinematic cuts that 
would be imposed in such a search. We apply those efficiencies to the inclusive tigg rate, but we treat 
the gluinos as final state particles, impose a minimum transverse momentum cut on the gluinos, and 
assume that their decay leads to an observable hadronic jet with vertex tag from the daughter bottom 
quark. Kinematic cuts used at the Tevatron (LHC) are as follows: 

p T (j) > 15(20)GeV, \ V (j)\< 3.0(4.0), 

pr{b) > 20(20)GeV, \n(b)\ < 2.0(2.5) , 
p T (l) > 15(15)GeV, < 2.0(2.5) , 

i> T > 30(30)GeV, AR mn > 0.4(0.4) . (1) 

where m, n are leptons, bottom quarks, gluinos or light jets. 

Matrix elements were constructed with a LESS-modified version of madgraph [H. We used 



CTEQ5L parton distribution functions [10] with factorization scale fif = m t + rrijj/2, where rrijj is 
the invariant mass of the extra bottom quark or gluino pair. The renormalization scale was taken to 
be the same, [i r = fif. We do not consider any additional contribution from tibibi production, as first 
the rate is much lower than for gluinos, and second as this introduces additional model dependence: 
whether the bottom squarks are long-lived or not; and if not, whether they have sufficient mass to 
decay into bottom quarks. 

In the LESS scenario, we take into account the altered running of the QCD coupling, a s , which 
occurs due to the presence of the light gluino and bottom squark contributions to the beta function. 
This causes a s to be considerably larger at the top quark mass scale. Since the cross sections are 
proportional to a A s) this effect increases the signal considerably. In calculating the signal cross section, 
we must also consider the effect of enhanced a s on the SM ttbb rate - an increase of 30%. We fix 
the value of the coupling to be a s (rrib) = 0.205, and use two-loop running with the bottom squark 
mass set to the bottom quark mass, m~ h = ||. We run the coupling from low-energy data, 
because if light gluinos and squarks exist, then the extraction of a s from low energy data would be 
unaffected by the LESS particle content. The value of a s obtained at the Z mass scale is 0.127-0.128 
for m g = 15 — 12 GeV. Given the uncertainty on a s (m{,), this result is within experimental errors, as 
discussed in Ref. ||. 

We calculate the ttbb cross section at the Tevatron to be 4.0 fb for prib) > 20 GeV (applied to 
the additional bottom quarks only; no top decays) []. We find a 25% efficiency for the kinematic cuts 
for both the semileptonic (brancyhing ratio (BR) = 29%) and all-hadronic (BR = 46%) decay modes 
of the top quarks. We do not consider the all-leptonic channel (BR = 4.7%) as the rate is much less 
than one expected event for reasonable Run II luminosity. The LESS cross section for tigg production 



1 Ref [11] imposed the same prib) cut but did not impose a cut on the rapidity of the b quarks. 
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varies from 11.2 fb at rrig = 12 GeV to 8.1 fb at 16 GeV, as shown in Fig. [T]. The tibb rate is 5.8 fb with 
LESS a s running. Assuming that the all-hadronic mode could be used (we don't assume top quark 
reconstruction is necessary), which is somewhat optimistic, and with a bottom quark vertex tagging 
efficiency of ej, = 50% and demanding that all four bottom quarks are tagged, then with 30 fb' 1 of 
integrated luminosity, we estimate between 4.5 and 3.5 signal events on a SM background Q of 1.4 
events. Using Poisson statistics, this corresponds to a 2.6 to 2.0 sigma effect. Alternatively, one may 
interpret this as the Tevatron having some capability to place 95% c.l. limits on this scenario. 

We note that this analysis is different from the planned search for tiH, which has the same final 
state signature, except for lack of a mass peak in the extra bottom quark pair spectrum. The fact that 
this final state arising from light gluinos produces essentially identical kinematic distributions to the 
SM case makes our proposed channel search more difficult, but is mitigated by the much larger overall 
rate from light gluinos that from a Higgs boson. We speculate that this search could be improved at 



the Tevatron by requiring only three vertex tags, as in the analysis of Ref. [|TTJ. This would increase 
the total sample by about a factor of five, but also approximately double the SM background by fake 
tags from tigg events. Nevertheless, 22 signal events on a background of about 15 events is an m 4.5a 
effect. We feel more thorough investigation along these lines is warranted. It may also be useful to 
examine additional production channels at the Tevatron, such as Zgg, which will have larger cross 
sections but other complications in their analysis. 

The situation is much better at the LHC We calculate the SM tibb total cross section {jprip) > 
20 GeV, no top quark decays) to be 1.9 pb, a rate copious enough to allow one to examine the cleaner 
all-leptonic top quark decay channel as well as the semi-leptonic channel. We find cuts efficiencies of 
30% and 20%, respectively. The tigg cross section varies from 9.0 pb for 12 GeV gluinos to 7.2 pb 
for 16 GeV gluinos, as shown in Fig. |I|. The tibb rate is 2.8 pb with LESS a s running. Again using 
€b = 50% and demanding four tags, we estimate that each experiment would observe from 15 to 12 
signal events in the all-leptonic channel alone, with only 2 fb -1 of data. Against the SM background 
of 3.3 events, this would yield a 5a effect over the entire gluino mass range considered. In the semi- 
leptonic channel for the same amount of data, we estimate 61 to 54 signal events on a background of 
18 events, potentially resulting in better than 12<r observation. 

One caveat is that of the long-lived bottom squark scenario: if the daughter bottom quark jets 
coming from the low-pr portion of the gluino spectrum do not have sufficient energy to be identified 
as jets, this analysis could suffer. To compensate we also investigated the case where all b partons and 
gluinos were instead required to have pr > 50 GeV. The rates fall by somewhat more than 50%, but 
the all-leptonic decay channel would need only 10 fb -1 and the semi-leptonic channel only 3 fb" 1 each 
to reflect a 5a observation of light gluinos. Thus, we feel this this analysis can easily be made model 
independent at the LHC. 

A final point to consider is the scale uncertainty of the cross sections. We estimate this by first 
calculating the cross sections for \i^ r = 2/i/ ir , yU/, r /2 and find about +75%/ — 45% variation. We also 
find a 75% enhancement if we use fif as originally stated, but apply two factors of a s (fi r = /i/) and 
two factors a s (rrijj). At LHC energies this results in an increase of about 75% in the rate. Since our 
signal cross sections are typically a factor four larger than the SM background, we are safe in assuming 
that theoretical uncertainties on the cross section cannot be a limiting factor in this analysis. 



2 Calculated with SM running of the coupling. 
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Figure 1. tigg cross section (solid lines) as a function of the gluino mass, m g , at the Fermilab Tevatron 
(left) and CERN LHC (right) for p T (b) > 20 GeV. Also shown are the SM tibb rates (dashed lines), and tibb 
rate with LESS-enhanced running of a s (dotted lines). 



We have presented an alternative production channel for light gluinos in a particular LESS model. 
It has the advantages of extremely distinctive final state signature (four bottom quarks and two W 
bosons), scale uncertainties much smaller than the increase in rate due to non-SM particles, and the 
only model dependence is the gluino mass. The observable rate at the Tevatron, after kinematical 
cuts and approximate efficiencies, is unfortunately useful only to place probably 95% c.l. limits on the 
LESS scenario, and even then only with large integrated luminosity. However, the LHC can make a 
5o" observation of light gluinos with only a few months of running at planned luminosity. 



ACKNOWLEDGMENTS 

We want to thank Alex Kagan for useful discussions. Fermilab is operated by URA under DOE 
contract No. DE-AC02-76CH03000. 



4 



Bibliography 



[1] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 71, 500 (1993); 79, 572 (1997); 75, 1451 
(1995); B. Abbott et al. [DO Collaboration], Phys. Lett. B 487, 264 (2000); B. Abbott et al. [DO 
Collaboration], Phys. Rev. Lett. 85, 5068 (2000); 

with certain changes in the theoretical input, the excess can be reduced, see M. Cacciari and 
P. Nason, |arXiv:hep-ph/ 02040251 
[2] P. Nason et al., Proceedings of the 1999 CERN Workshop on Standard Model Physics (and more) 
and the LHC, ed. G. Altarelli and M. L. Mangano (CERN, Geneva, 2000), p. 231; S. Frixone, 
M. L. Mangano, P. Nason and G. Ridolfi, Heavy Flavors II, ed. A. J. Buras and M. Linder (World 
Scientific, Singapore, 1997). 



[3] 
[4] 
[5] 
[6] 
[7 
[8 
[9 
[10 

[11 



E. L. Berger et al, Phys. Rev. Lett. 86, 4231 (2001). 

H. P. Nilles, Phys. Rept. 110, 1 (1984); H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985). 

See, however, J. j. Cao, Z. h. Xiong and J. M. Yang, |arXiv:hep-ph/ 01111 44j . 

M. Carena, S. Heinemeyer, C. E. Wagner and G. Weiglein, Phys. Rev. Lett. 86, 4463 (2001). 

E. L. Berger and L. Clavelli, Phys. Lett. B 512, 115 (2001) ||arXiv:hep-ph/0105147 l. 

T. Becher, S. Braig, M. Neubert and A. Kagan, arXiv: hep-ph / 1 1 2 1 29 . 

T. Stelzer and W. F. Long, Comp. Phys. Comm. 81, 357 (1994). 

H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12, 375 (2000). 

J. Goldstein et al, Phys. Rev. Lett. 86, 1694 (2001). 



5 



